HU-EP-08/59, ITEP-L AT/2008-22 



Infrared behavior and Gribov ambiguity 
in SU{2) lattice gauge theory 

V. G. Bornyakov 
Institute for High Energy Physics , 142281, Protvino, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 



V. K. Mitrjushkin 
Joint Institute for Nuclear Research, 141980 Dubna, Russia 
and Institute of Theoretical and Experimental Physics, 117259 Moscow, Russia 



00 

o 

O 
(N 

CJ 

Q 

in 



oo 

o 



M. Miiller-Preussker 
Humboldt Universitdt zu Berlin, Institut fiir Physik, Newton-Str. 

(Dated: December 15, 2008) 



15, 12489 Berlin, Germany 



For SU (2) lattice gauge theory we study numerically the infrared behavior of the Landau gauge 
ghost and gluon propagators with the special accent on the Gribov copy dependence. Applying a 
very efficient gauge fixing procedure and generating up to 80 gauge copies we find that the Gribov 
copy effect for both propagators is essential in the infrared. In particular, our best copy dressing 
function of the ghost propagator approaches a plateau in the infrared, while for the random first 
copy it still grows. Our best copy zero-momentum gluon propagator shows a tendency to decrease 
with growing lattice size which excludes singular solutions. Our results look compatible with the 
so-called decoupling solution with a non-singular gluon propagator. However, we do not yet consider 
the Gribov copy problem to be finally resolved. 
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I. INTRODUCTION 

The lattice study of the gluon and ghost propaga- 
tors in Landau gauge has a long history. It was started 
by Mandula and Ogilvie [l| and continued in many pa- 
pers (for review see, e.g., and references therein). 
One of the main goals of such studies was to clar- 
ify the infrared (IR) asymptotics of the propagators 
and of the running coupling which can be determined 
through these propagators. The hope always was that 
ab-initio lattice results will give support to or discrimi- 
nate between various theoretical predictions for the IR 
behavior obtained with continuum methods, in par- 
ticular within the Dyson-Schwinger approach. One 
prediction was definitely overturned: lattice results 
showed that the gluon propagator is not divergent in 
the infrared. 

At the same time it has been found that the lat- 
tice approach has its own difficulties when applied to 
such studies. One of them is that to reach small mo- 
menta necessary to study the IR limit one has to go 
to huge lattices which makes the numerical simula- 
tions formidable. Another, less apparent, but not less 
difficult problem is the problem of Gribov copies. Al- 
though for many years it was believed that the effect 
of Gribov copies on both gluon and ghost propagators 
was weak and could be considered just as a noise in the 
scaling region 0,13] it has been found first for the ghost 
propagator and quite recently for the gluon prop- 
agator that these effects are in fact quite strong. 
The presence of these effects makes the task of lattice 



computations of the field propagators in the IR region 
even more difficult. 

These difficulties of the lattice approach has made it 
impossible so far to obtain results which could confirm 
or disprove the existent confinement scenarios pro- 
posed by Gribov 0] and Zwanziger Q on one hand 
and by Kugo-Ojima 0] on the other. The Gribov- 
Zwanziger scenario predicts that the gluon propaga- 
tor is IR vanishing, while the Kugo-Ojima criterion of 
confinement predicts the ghost dressing function to be 
IR divergent. 

In recent years the interest to the lattice results 
for the field propagators in the IR region has been 
revived. This interest was stimulated also by the 
practical progress achieved over the years within 
the Dyson-Schwinger (DS) approach as pursued by 
Alkofer, von Smekal and others (for an intermedi- 
ate review see [l^), and more recently with the help 
of functional renormalization group (FRG) equations 
[nl, Infrared QCD has been also investigated us- 
ing the stochastic quantization method [ij, as 
well as with effective actions [H, [l^ . 

In this paper we continue our lattice study of the 
influence of Gribov copies on the (minimal) Landau 
gauge SU{2) gluon and ghost propagators in the IR 
region by applying global Z{2) flip transformations 
in combination with an effective optimization algo- 
rithm, the so-called simulated annealing (SA). The flip 
transformation was introduced in |f| . Its infiuence on 
the gluon propagator was thoroughly studied lateron 
fVT]. The Z{2) flips - equivalent to non-periodic Z{2) 
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gauge transformations - were shown to cause rather 
strong effects in the IR behavior of the gluon prop- 
agator. In [l3| high statistics computations of the 
gluon propagator were made for lattice sizes varying 
from 1.8 fm to 6.5 fm at one fixed bare lattice cou- 
pling (3 = A/ — 2.20. The latter was chosen in order 
to reach reasonably large physical volumes (and thus 
small momenta) on comparatively moderate lattice 
sizes up to 32''. It turned out that due to better gauge 
fixing finite-volume effects, usually strong at minimal 
momenta, became largely suppressed. Furthermore, it 
has been observed that at momenta p ~ 270 Mev the 
gluon propagator seems to have a turning point leav- 
ing open the possibility for a vanishing gluon propa- 
gator in the IR limit ^ 0. Here we continue this 
investigation enlarging the lattice up to 40^* at the 
same /3-value corresponding to a volume of (8.4fm)'' 
and extending the studies to the ghost propagator, 
too. We systematically search for Gribov copies by 
combining all 2* — 16 Z{2) Polyakov loop sectors for 
all Euclidean directions into one gauge orbit. 

The main motivation for this computation is trig- 
gered by the puzzle posed by the above mentioned 
continuum approaches. Different kinds of solutions 
with a quite different IR behavior of the gluon and 
ghost propagators have been reported by different 
groups. The power-like solution with relation between 
gluon Ku and ghost kg > exponents kd = — 2kg 
was called recently a scaling solution [l^. This solu- 
tion [H, [H, [H, 0, [m allows the gluon propagator 
to vanish and the ghost dressing function diverge in 
the IR limit in one-to-one correspondence with both 
the Gribov-Zwanzi ger scenario [7|, |22] and the Kugo- 
Ojima criterion (ol. l23l| for confinement. On the con- 
trary, the so-called decoupling solutions [l^. [2^. [25l. [26t 
provide an IR finite or weakly divergent gluon propa- 
gator and a finite ghost dressing function leading to a 
running coupling vanishing in the infrared. For recent 
discussions of the present status of research see, e.g., 
[ist and further references therein'. The lattice ap- 
proach based on the first-principle path integral quan- 
tization should be able to resolve the issue. 

Results for SU{2) [13 as well as for SU{i) 28] ob- 
tained on very large lattices and by employing purely 
periodic gauge transformations seem to be in conflict 
with the scaling solution and are compatible with the 
decoupling solution. That this might be not in conflict 
with the (appropriately modified) Gribov-Zwanziger 
scenario has been recently pointed out in (29| . 

Here, by enlarging the gauge orbits with non- 
periodic Z{2) flip gauge transformations and employ- 
ing the SA algorithm we shall come closer to the global 



extremum of the Landau gauge functional, i.e. closer 
to the fundamental modular region. We find the Gri- 
bov copy dependence to be very strong. Still our re- 
sults look rather as an argument in favor of the decou- 
pling solution with a non-singular gluon propagator. 
However, we do not yet consider the problem of Gri- 
bov copies and, correspondingly, the infrared asymp- 
totics of the gluon propagator to be finally resolved. 

In Section [H] we introduce the observables to be 
computed. In Section HTTl some details of the gauge fix- 
ing method and of the simulation are given, whereas 
in Section lTVl we present our results. Before coming to 
the conclusions in Section IVTl we will discuss the de- 
pendence of our results on the number of gauge copies 
in Section Ivl 



II. GLUON AND GHOST PROPAGATORS: 
THE DEFINITIONS 

For the Monte Carlo generation of ensembles of non- 
gauge-fixed gauge field configurations we use the stan- 
dard Wilson action, which for the case of an SU{2) 
gauge group is written 



1 - i Tr 
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Uxfj.Ux + fj,-iyUl_^^.^UlJj 



(1) 



Here gQ is a bare coupling constant and U^^ G SU(2) 
are the link variables. The latter transform as follows 
under gauge transformations gx 

C/.M ^ C^.'m = 9lUx,,9x+^. , gx e SU{2) . (2) 

The standard definition [ij of the dimensionless lattice 
gauge vector potential Ax+ii/2,iJ, is 



jAx 



The reader should keep in mind that the definition 
is not unique which can have an essential influence 
on the propagator results in the IR region, where the 
continuum limit is hard to control. 

In lattice gauge theory the usual choice of the Lan- 
dau gauge condition is [l|] 



idA)x = J2 ('4^+/'/2;P - A-/./2;m) = , (4) 

which is equivalent to finding an extremum of the 
gauge functional 



^ It is worthwhile to note that the DS approach introduced 
originally as a method for resummation of the perturbative 
series is not sensible to different gauge copies. 
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where V — is the lattice volume, with respect to 
gauge transformations gx ■ After replacing U ^ 
at the extremum the gauge condition @ is satisfied. 
The manifold consisting of Gribov copies providing 
local maxima of the functional (O and a semi-positive 
Faddeev-Popov operator (see below) is called Gribov 
region Q, while that of the global maxima is called 
the fundamental modular domain A C fi. Our gauge 
fixing procedure is aimed to approach this domain. 

The gluon propagator D and its dressing function 
Z are then defined (for p ^ 0) by 



Zip) = D{p) , (6) 

where A{k) represents the Fourier transform of the 
gauge potentials defined by Eq. ([3]) after having fixed 
the gauge. The momentum p is given by p^ = 
{2 / a) sin {nk^/L), k^ e {-L/2,L/2]. For p / 0, 
one gets from Eq. ^ 



a— 1 1 

whereas at p — the "zero momentum propagator" 
D{0) is defined as 

a— 1 /^— 1 

The lattice expression for the Landau gauge 
Faddeev-Popov operator Af"'' = -d^,Df (where Df 
denotes the covariant derivative in the adjoint repre- 
sentation) for SU{2) is given by 

Mfy[U]=Y^ {{Sf, + Sf_^,,) 5x;y 

~ {,^xn ~ 5y-x+fi (9) 

where 



qab cab rp jj a ah ahc ac 

'-'xfi — " 2 ' ~ ~2 x+fi/2;tj. ■ 

(10) 

From the expression (fTO|) it follows that a trivial zero 
eigenvalue is always present, such that at the Gri- 
bov horizon dQ the first non-trivial zero-eigenvalue 



appears. Thus, if the Landau gauge is properly imple- 
mented, M[U] is a symmetric and semi-positive ma- 
trix. 

The ghost propagator G°-^{x,y) is defined as [sl-lsojl 

G^\x,y) = cJ"" G{x - y) ^ ^( )l\ [U] ) ■ 

(11) 

Note that the ghost propagator becomes translational 
invariant [i.e., dependent only on x — y) and diagonal 
in color space only in the result of averaging over the 
ensemble of gauge-fixed representants of the original 
Monte Carlo gauge configurations. The ghost propa- 
gator G{p) in momentum space and its dressing func- 
tion J{p) can be written as 

X, y, a 

J(p) = G{p) p' , (12) 

where the coefficient is taken for a full normal- 
ization, including the indicated color average over 
a = 1,2,3. We mentioned above that M[U] is sym- 
metric and semi-positive. In particular, it is positive- 
definite in the subspace orthogonal to constant vec- 
tors. The latter are zero modes of M[[/]. Therefore, it 
can be inverted by using a conjugate-gradient method, 
provided that both the source il}°'{y) and the initial 
guess of the solution are orthogonal to zero rnodes. 
As the source we adopted the one proposed in [4] 

^l;''{y) ^ 5'"'e^'''P-y (0,0, 0,0), (13) 

for which the condition ip°'{y) = is automati- 
cally imposed. Choosing the source in this way allows 
to save computer time since, instead of the summa- 
tion over X and y in Eq. (|12p . only the scalar product 
of M^^ip with the source ip itself has to be evaluated. 
In general, the gauge fixed configurations can be used 
in a more efficient way if the inversion of M is done 
on sources for c — 1,2,3 such that the (adjoint) color 
averaging, formally required in Eq. p2p . will be ex- 
plicitely performed. 

III. SIMULATION DETAILS 

We restrict ourselves to Monte Carlo (MC) simu- 
lations at f3 = 4/^0 — 2.20 and use lattice field con- 
figurations for which the gluon propagator has been 
already computed in [l7l|. Here we add the computa- 
tion of the ghost propagator and new data obtained 
on the larger symmetric lattice with the linear size of 
L = 40. For the latter case we generated an ensemble 
of 430 independent Monte Carlo lattice field configu- 
rations. Consecutive configurations (considered as in- 
dependent) were separated by 100 sweeps, each sweep 
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TABLE I: Lattice sizes, statistics, number of gauge copies 
used throughout this paper. The second (third) column 
gives the number of configurations used to compute the 
gluon (ghost) propagators. 



being of one local heatbath update followed by L/2 
microcanonical updates. In TableUwe provide the full 
information about the field ensembles used through- 
out this paper. 

For gauge fixiiig we employ the Z{2) flip operation 
as discussed in [131 • It consists in flipping all link 
variables Ux^_i attached and orthogonal to a Sd plane 
by multiplying them with —1. Such global flips are 
equivalent to non-periodic gauge transformations and 
represent an exact symmetry of the pure gauge action 
considered here. The Polyakov loops in the direction 
of the chosen links and averaged over the 3ci plane 
obviously change their sign. Therefore, the flip oper- 
ations combine for each lattice field configuration the 
2"* distinct gauge orbits (or Polyakov loop sectors) of 
strictly periodic gauge transformations into one larger 
gauge orbit. 

The second ingredient of our gauge fixing procedure 
is the consequent use of the simulated annealing (SA) 
method, which has been found even computationally 
more efficient than the only use of standard overrelax- 
ation (OR). 

The SA algorithm generates a field of gauge trans- 
formations g{x) by MC iterations with a statisti- 
cal weight proportional to eyi\i{W Fu[g]/T) . The 
"temperature" T is a technical parameter which is 
gradually decreased in order to maximize the gauge 
functional Fu [g] . In the beginning, T has to be cho- 
sen sufficiently large in order to allow traversing the 
configuration space of g{x) fields in large steps. It has 
been checked that an initial value Tinit = 1.5 is high 
enough. After each quasi-equilibrium sweep, includ- 
ing both heatbath and also microcanonical updates, 
T has been decreased with equal step size until g{x) 
is uniquely captured in one basin of attraction. The 
criterion of success is that during the consecutively 
applied OR the violation of transversality decreases 
in a more or less monotonous manner for almost all 
applications of the compound algorithm. This condi- 
tion turned out reasonably satisfied for a final lower 
temperature value Tfinai — 0.01 [sj]. The number 
of temperature steps was chosen to be 1000 for the 



smaller lattice sizes and has been increased to 2000 
for the lattice size 40^ included here. The finalizing 
OR algorithm requires a number of iterations varying 
from 0(10^) to 0(10^). In what follows we shall call 
the combined algorithm employing SA (with finalizing 
OR) and Z{2) flips the 'FSA' algorithm. 

Some details of the gauge fixing procedure com- 
pared to our previous work p^j have been changed. 
For every configuration the Landau gauge was fixed 
Ncopy = 80 times (5 gauge copies for every flip- 
sector), each time starting from a random gauge 
transformation of the mother configuration, obtain- 
ing in this way Ncopy Landau-gauge fixed copies. In 



[17[ where smaller lattices were simulated Ncopy was 
smaller: 40 (for 32" lattice), 24 (24"* and 16*). Only on 
very small lattices 12"* and S'', where producing copies 
was substantially cheaper, we produced 80 copies. 

In order to reduce the computational effort in the 
finalizing OR sweeps on the 40* lattice we applied the 
following trick. We noticed that after a comparably 
small number of OR sweeps, definitely before the con- 
vergence criterion is reached, one can already decide 
which copy has a higher maximum of the gauge func- 
tional, i.e. one can stop the OR procedure already 
when the change in the functional becomes compara- 
bly small and further sweeps will not change the order 
of copies according to the value of the maximized func- 
tional. Note, that the functional was different from its 
final value only in the 8th digit and we used these val- 
ues in Fig. [71 After having selected the 'best copy' 
[he ) the OR gauge fixing for this copy has to be final- 
ized. To be precise, complete gauge fixing was made 
also on the randomly chosen 'first copy' (/c ), just for 
the purpose of comparison. 

For the finalizing OR we used the standard Los 
Alamos type overrelaxation with the parameter oj = 
1.7. For he and /c the iterations have been stopped 
when the following transversality condition was satis- 
fied: 



max 

x, a 



(14) 



We used the parameters eior = 10 (i.e. 10 for 
(9A)2). 



IV. RESULTS 

In this section we present the data for the gluon and 
ghost propagator. In Fig. [1] we show the new data for 
the gluon propagator D{p) in physical units obtained 
on the 40* lattice at /? = 2.20. We compare the be FSA 
result with the fe SA result (the latter without fiips). 
We clearly see the Gribov copy effect for the lowest 
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FIG. 1: The momentum dependence and Gribov copy sen- 
sitivity of the gluon propagator D{p) in the IR region on 
the 40* lattice. FiUed symbols correspond to the be en- 
semble, open symbols to the fc ensemble. 



> 
CD 



12 
11 
10 
9 
8 
7 
6 
5 
4 



D(0), FSA - 








U(UJ, UK ^ 


— O 1 












* 


# * 

















0.02 



0.04 0.06 
1/L, GeV 



0.08 



0.1 



FIG. 3: The dependence of D{0) on the lattice size. 
be FSA results are compared with fe OR results (without 
flips). 



7 

> 

o 5 

-a 4 

° 3 
2 
1 




* 8- 




V 12 




° 16* 




« 24* 




o 32* 




• 40* 





0.0 0.5 1.0 1.5 2.0 
Ipl; Gev 



2.5 



3.0 



FIG. 2: The momentum dependence of the gluon prop- 
agator D{p) on various lattice size, be results are shown 
throughout. 



accessible momenta moving the data points to lower 
values for better copies (with the larger gauge func- 
tional) . The different points at p ~ SOOMeV belong to 
different realizations of and seem to indicate some 
violation of the hypercubic symmetry. 

In Fig. [2] we present these new data together with 
the ones obtained on smaller lattice sizes always for 
the FSA be case. 

We see that the data are nicely consistent with each 
other and indicate a turnover to decreasing values to- 
wards vanishing momentum. A smooth extrapola- 
tion to D{0) becomes visible. But still there is no 
indication for a vanishing gluon propagator at zero 
momentum for increasing volume. This is demon- 
strated in Fig. [21 where we show the dependence of 



the zero- momentum propagator -D(O) as a function 
of the inverse linear lattice size 1/L. This behav- 
ior demonstrates a (slight) tendency to decrease, and 
looks hardly consistent with D{0) = limit. One 
could consider it rather like an argument in favor of 
the decoupling solution with a finite gluon propaga- 
tor in the infrared. However, one still cannot exclude 
that there are even more efficient gauge fixing meth- 
ods, superior to the one we use, which could make this 
decreasing more drastic. 

Using Ward-Slavnov- Taylor identities the authors 
of [24, 26, 32] came to the conclusion that the gluon 
propagator should be IR divergent, however, this di- 
vergence might be so weak that it could be hardly 
resolved on the lattice. We believe that our results 
for D{0) are in clear disagreement even with a weak 
divergence. 

Analogously to Fig. [T] in Fig. |4] wc show the ghost 
dressing function J{p) obtained on the 40^ lattice. 
There is a very clear Gribov copy effect changing J{p) 
even qualitativly. Whereas the fc SA results seem 
to support a weakly singular behavior, the be FSA 
data provide a plateau pointing to a finite IR value of 
the ghost dressing function, i.e. a tree-level behavior 
of the ghost propagator. Our data indicate that the 
plateau starts at p < 200 MeV. 

In Fig. O the ghost dressing function is shown for 
lattice sizes from 16^ to 40'*. We show always be FSA 
results, except for 24^, where we compare also with 
fc data obtained with the conventional OR algorithm. 
The latter show an even stronger IR singular behavior 
than those data obtained with the fe SA algorithm. 

There is a clear weakening of the singularity visi- 
ble additionally to a finite-size effect which seems to 
lead to an IR plateau behavior. Such a plateau would 
be consistent with the different decoupling solutions 
and in contradiction with the Kugo-Ojima confine- 
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such that an approach to an IR fixed point as expected 
from the scaling DS and FRG sohition seems to be 
excluded. 
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ment criterion [2^ . 

In Fig. [6] for the be FSA results obtained on lat- 
tice sizes from up to 40^ we draw the behavior of 
the running coupling related to the ghost-ghost-gluon 
vertex 



^sip) - §^J'{P) Z{p) 



(15) 



under the assumption that the vertex function is con- 
stant as seen in perturbation theo ry 1 331 and approx- 
imately also in lattice simulations [34l . |35| . 

The decrease towards = is obvious. With the 
improved gauge fixing the effect is even strengthened, 



V. DISCUSSION: THE QUALITY OF THE 
GAUGE FIXING PROCEDURE 

In most of our simulations we have generated up to 
Ncopy — 80 gauge copies for every thermalized con- 
figuration (up to 5 gauge copies for every flip-sector). 
A very reasonable question is whether our results will 
change if we will further increase the number of gauge 
copies Ncopy In Fig. [7] we show the dependence of the 
average be functional {Fi,c){kcopy) ■ 



1 " 

(-^bc)(^copy) — ^ ^ ^bc ij^copy) ; ^copy ~ 1 1 



.N, 



copy 5 

(16) 

where for every configuration Fi,c{kcopy) is the 'best' 
(i.e., maximal) value of the functional F found after 
employing kcopy copies, and n denotes the number of 
configurations. 

One can see that this average still keeps a tendency 
to grow, which could mean that one should take even 
more copies to reach the global maxima. To under- 
stand it better we generated 25 configurations on the 
40^ lattice with Ncopy = 320 (i.e., 20 gauge copies per 
sector). 

In Fig. [5] we show the difference 

AF{kcopy) = Fbcikcopy) " Fbc{80) (17) 

at kcopy — 160, 240 and 320 for these configurations. 

For the majority of configurations this difference is 
rather small. However, for about 20% of configura- 
tions the difference is of the order of 10~^, which could 
mean still a rather strong influence on the values of 
propagators. 
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To demonstrate that the change in the functional of 
the order of 10~^ indeed might give rise to substantial 
change in the propagators we plot in Fig. [9] the gluon 
propagator computed on 32^ lattice at two momenta, 
p = and p = Pmin as a function of the difference 
Fbc{32) - Fbcikcopy) for kcopy = 1,2,..., 18. One can 
see that the change in the functional Ff,c in the 5th 
digit brings quite substantial change in the propagator 
at both momenta. 

These observations give an idea that there might be 
another even more efHcient gauge fixing method which 
will be more successful in the search of the global max- 
imum of the functional F (as FSA is superior with 
respect to standard OR). 

Therefore, we cannot draw a final conclusion before 
spending much more efforts into optimization of the 
gauge fixing procedure. However, we expect that both 
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the gluon propagator as well as the ghost propagator 
will become further suppressed in the infrared, when 
approaching the fundamental modular region. 



VI. CONCLUSIONS 

In this work we studied numerically the dependence 
of the Landau gauge gluon and ghost propagators, as 
well as of the running coupling constant, in pure gauge 
SU (2) lattice theory in the infrared region. The spe- 
cial accent has been made on the study of the depen- 
dence of these 'observables' on the choice of Gribov 
copies. 

The simulations have been performed using the 
standard Wilson action at (3 = 2.20 for linear lat- 
tice sizes up to i = 40. For gauge fixing gauge orbits 
enlarged by Z(2) flip operations were considered with 
up to 5 gauge copies in every flip-sector (in total, up 
to 80 gauge copies). For 25 thermalized configurations 
we produced 20 copies per sector (in total, 320 gauge 
copies for every configuration). The maximization of 
the gauge functional was achieved by the simulated 
annealing method always combined with consecutive 
overrelaxation ('FSA' algorithm). 

Our findings can be summarized as follows. 

1) For the gluon propagator our new data for the 
40^ lattice agree with data on the smaller lattices (up 
to 32^*). We confirm our conclusion [l3| about the 
appearance of the local maximum at a non-zero value 
of the momentum (this local maximum was absent 
for lattice sizes < 24"*). 

The zero-momentum gluon propagator D{0) has a 
tendency to decrease with growing lattice size L. This 
observation is in clear contradiction with the infrared 
divergent gluon propagator obtained on the basis of 
Ward-Slavnov- Taylor identities. 

At the time being, this behavior looks hardly con- 
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sistent with a £'(0) ~ limit at infinite L, and could 
be considered rather like an argument in favor of the 
decoupling solution with a non-singular gluon propa- 
gator. However, we do not yet consider the problem 
of the infrared asymptotics of the gluon propagator as 
a finally resolved (see below). 

2) We calculated the ghost propagator for lattices 
up to 40^^. Our be dressing function J{p) of the ghost 
propagator demonstrates the approach to a plateau in 
the infrared, while the fc dressing function still grow s 
(as it was in earlier calculations; see, e.g., [sl. [28l. 136} ) . 

This is a first clear indication of the lack of the IR- 
enhancement of the ghost propagator. This plateau 
behavior is in a clear contradiction with the Kugo- 
Ogima confinement criterion. The fate of this con- 
finement criterion still needs a further clarification. 

3) We have found that the effect of Gribov copies 
for both the propagators and in the consequence for 
the running coupling is essential in the infrared range 
p < 1 GeV. Therefore, the quality of the gauge fixing 
procedure in the study of gauge dependent observables 
remains important. 

Indeed, the FSA method provides systematically 
higher values of the functional Fu{g) as compared 
to the standard OR procedure for the same thermal- 
ized configurations. This means that in practice OR 
needs many more random copies to explore (corre- 
spondingly, much more CPU time to spend) to find 
larger values of Fu{g) as compared to FSA. This effect 
becomes stronger with increasing the volume. How- 
ever, we cannot say that we have reached the funda- 
mental modular region when fixing the Landau gauge 
on larger lattices. One cannot exclude that there is 
another method superior to our FSA algorithm. We 
believe that the Gribov problem deserves even more 



thorough studies. 

Maybe, there are alternative ways to resolve the 
problem of the IR asymptotics of the propagators. We 
have to be aware that the lattice method as it is nor- 
mally used has still some uncertainties. First of all 
the continuum limit in the infrared range is hard to 
control and it depends on the proper choice of the 
gauge potential Ax^- Moreover, the infrared limit 
is sensitive to the boundary conditions, which nor- 
mally are taken to be periodic. Incomplete gauge fix- 
ing in combination with these choices seems unavoid- 
ably to lead to zero-momentum modes not sufficiently 
suppressed even in the thermodynamic limit. That 
the presence of zero-momentum modes can spoil the 
behavior of gauge-variant propagators is well-known 
from the example of 4d compact C/1) lattice gauge 
theory [13, H Isl- Whether a BRST conformal lat- 
tice reformulation will solve the issue as proposed in 
[iol . im remains to be seen. 
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